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Stream salamanders in the family Plethodontidae constitute a large 
biomass in and near headwater streams in the eastern United States and are 
promising indicators of stream ecosystem health. Many studies of stream 
salamanders have relied on population indices based on counts rather than 
population estimates based on techniques such as capture-recapture and 
removal. Application of estimation procedures allows the calculation of 
detection probabilities (the proportion of total animals present that are 
detected during a survey) and their associated sampling error, and may be 
essential for determining salamander population sizes and trends. In 1999, 
we conducted capture-recapture and removal population estimation 
methods for Desmognathus salamanders at six streams in Shenandoah 
National Park, Virginia, USA. Removal sampling appeared more efficient 
and detection probabilities from removal data were higher than those from 
capture-recapture. During 2001-2004, we used removal estimation at eight 
streams in the park to assess the usefulness of this technique for long-term 
monitoring of stream salamanders. Removal detection probabilities ranged 
from 0.39 to 0.96 for Desmognathus, 0.27 to 0.89 for Eurycea and 0.27 
to 0.75 for northern spring (Gyrinophilus porphyriticus) and northern red 
(Pseudotriton ruber) salamanders across stream transects. Detection 
probabili 
but did differ among streams for Desmognathus. Population estimates of 
Desmognathus decreased between 2001-2002 and 2003-2004 which may 
be related to changes in stream flow conditions. Removal-based procedures 
may be a feasible approach for population estimation of salamanders, but 
field methods should be designed to meet the assumptions of the sampling 
procedures. New approaches to estimating stream salamander popul 
are discussed. 
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INTRODUCTION 


Stream salamanders including Desmognathus, Eurycea, Gyrinophilus and Pseudotruon 
species in the family Plethodontidae (lungless salamanders) are found in and near seeps and 
streams in the eastern United States. These salamanders play an important role in nutrient 
cycling and energy flow and can be top vertebrate predators in fishless headwater streams 
(Davic, 2002). In Appalachian old growth forests, plethodontid salamanders are extremely 
abundant (up to 18,486 individuals/ha) and constitute a large biomass (16.53 kg/ha) exceeding 
that of birds (PETRANKA & MURRAY, 2001) Stream salamander larvae develop in seeps and 
streams, After transformation, juveniles and adults spend part of their hfe in the leaf litter, 
rocky substrate and banks along streams, foraging on the surface on wet or humid nights, and 
hiding beneath rocks, logs, leaves, moss, bark, and in burrows during the day (PETRANKA, 
1998) 


Because stream salamanders may serve as indicators of stream ecosystem health (CORN 
& Bury, 1989; PeTRANKA ct al., 1993. WELSH & OLLIVIER, 1998; SOUTHERLAND et al., 2004), 
identifying reliable survey methods and population estimation techniques for this group are 
important. Many studies of stream salamanders have used population indices such as raw 
counts and densities (WELSH & OLLIVIER, 1998, Barr & Bassitt, 2002), although studies 
using population estimates based on capture-recupture and removal sampling have also been 
conducted (BRUCE, 1995; Nut 15 & KAPLAN, 1998, PETRANRA & Murray, 2001). Population 
estimates tbat include detection probabilities (the proportion of total animals present that аге 
detected during a survey) and their associated samphng error may be essential in determining 
population sizes and trends. 


Detection probabilities ) for salamanders can vary spatially, temporally, and by species 
and age class (JUNG et al., 2000; SALvipio, 2001; Bau rv et al , 2004) If js differ among study 
sites or over time, population indices will not be comparable unless differences in detection are 
estimated during sampling Also, magnitude of detection probability can influence precision 
and bias of estimated population sizes; if p's or recapture rates of animals are low, standard 
errors of estimated population sizes become large, leading to unrchable or biased estimates. 


The purpose of this investigation was two-fold. A study conducted in 1999 was designed 
to compare capture-recapture and removal estimation methods for stream salamanders and 
determine which method was more efficient and provided higher detection probabilities, A 
study conducted from 2001-2004 was meant to assess the usefulness of the preferred tech- 
nique (removal estimation} for long-term monitoring of stream salamanders in Shenandoah 
National Park, Virginia, USA. 


MATERIALS AND METHODS, 


For all stream salamander surveys in Shenandoah National Park. we recorded salaman- 
der species, age class (those with gills were recorded as larvae and those without were classified 
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as adults), and snout-vent length (SVL) and total length (in mm) Because larvae of some 
species were occasionally not distinguished in the field and/or because of low sample sizes, we 
combined data for northern dusky (Desmognathus fuscus) and seal (D. monticola) salaman- 
ders into Desmognathus and combined data for northern spring (Gyrinophilus porphvriticus) 
and northern red {Pseudotriton ruber) salamanders into Gyrinophilusi Pseudotriton for analy- 
ses. The northern two-lined salamander (£urycea brsfineata) is found in the northern part of 
the park, whereas the southern two-Iied salamander (E curigera) is found in the southern 
part (GutrkA & SATLER, 1990). Herein, we treat the two species, which are indistinguishable 
morphologically, as one when streams throughout the park are considered. Stream salaman- 
der species in the park not detected in our surveys included the long-tailed (Е longicauda) and 
three-lined (E. guttolimeata) salamanders (Wrrr, 1993). 


In 1999, we compared capture-recapture and removal methods to estimate stream 
salamander populations at six streams in the park. Three of the streams were first order 
(Jeremy’s Run, Land's Run, Piney River) and the other three were second order. Two of four 
observers turned over the top layer of rocks and other objects greater than 6.4 cm maximum 
width or length within two 50 x 1 m transects (one on each side of the stream channel, 100 m? 
total area) at each stream, sampling the terrestrial habitat immediately adjacent to the wetted 
stream channel Larvae were captured in the hy porheic zone and comprised 9 + 2,8". (mean 
+ s) of Desmognathus observations, 35 + 8,7° of E bnlimneata observations, and 68 + 
11.6 % of Gyrinophilus! Pseudotriton observations across all surveys. 


Approximately weekly from 9 July to 13 August 1999, we captured salamanders during 
the day by hand and dipnet and batch marked larvae and adults of or above 25 mm SVL using 
visible implant fluorescent elastomer (VIE, Northwest Marine Technologies, Inc ). a biocom- 
patible latex-based dye injected just under the skin. We used one of three VIE colors (orange. 
red. green) at one of four positions just behind the forelimbs or in front of the hindlimbs and 
checked marks under a blanket using an ultraviolet light (JUNG et al. 2000) Other studies 
have shown that VIE marks are more permanent and cause salamanders less harm than 
toe-clipping (Davis & Ovaska, 2001; Maroro, 2001) Identities and measurements of 
unmarked, marked and recaptured salamanders were recorded at cach survey The percent of 
salamanders that escaped per stream across capture-recapture surveys averaged 36 2.2"; 
for Desmognathus, 38 + 55 » for E. Het. and 43 + 10.3 , for Giimaphilust 
Pseudotriton, 


After completing 5-6 capture recapture surveys (five at Keyser Run, Pass Run Tributary, 
Piney River), we conducted temporary removal sampling of stream salamanders from 23 to 
28 August 1999 at the same transects. Three passes, at least two hours apart, were made each 
day for two consecutive days for a total of six removal passes. We tallied the number of larvae 
and adults of each species removed at cach pass and kept species and age classes in separate 
buckets with stream water positioned in the stream in the shade All salamanders were 
released back to transects after the final pass. 


During June and July of 2001-2004, we used temporary removal sampling as above 
though with 2 3 passes conducted one after the other during the day at 1-2 transects at each of 
eight streams in the park Five of the streams were first order and three were second order or 
higher Surveys were also conducted at a ninth stream. Staunton River, but no salamanders 
were detected there, presumably due to residual elfects of a large flood along this river m 1995 
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For these surveys, we used a modified transect design, transects were 15 m long and 2 m wide 
and located on only one side of the stream, spanning I m along the stream bank and 1 m in the 
stream channel. This design allowed for capture of significantly more larval £urycea(F = 21.9, 
df = 1, P < 0001) and GyrmoplulusiPseudotriton (F = 4 4, df = 1, Р O 043) salamanders 
compared to the 1999 stream transects, larvae comprised 5 + 22° of Desmognathus 
observations, 75 = 4.1% of Ewrvcea observations, and 89 + 4.2 % of Gyrinoplulus/ 
Pseudotriton observations Our summer surveys mostly missed Desmognathus larvae, which 
typically transform by June-July (PETRANKA, 1998). 


For capture-recapture data from 1999. we estimated population sizes (N) and detection 
probabilities (5) and their standard errors (s,) using maximum likelihood and Bayesian 
estimators, These models assume a closed population and fit a series of models that differ in 
their assumptions about variation in 5 during sampling (Orts et al, 1978, Rexstap & 
BURNHAM, 1991). The test for population closure in program CAPTURE showed that all 
populations at cach stream were closed (all P > 0 05), so open population models were not 
used (Rrastap & Викмнам, 1991) At least one assumption of capture-recapture models, 
that all animals captured are marked, was not met because we estimated only a subset of the 
populations, i.e., individuals of or above 25 mm SVL The Bayesian estimator (Gazry & 
SrALL Y, 1986) assumes prior distributions for М and p and estimators are derived based on the 
posterior distribution, We used a umform (0,1) prior distribution on p and a diffuse negative 
binomial prior distribution on № (GrogGE & Кові кт, 1992) We fit the Bayesian model using 
a Markov chain Monte Carlo technique known as Gibbs sampling. in which the posterior 
distributions are estimated by simulation Bayesian estimators do not rely on the asymptotic 
properties of maximum likebhood estimators, and hence are preferred for small sample sizes 
(Gazey & STALEY, 1986) 


For removal data from 1999 and 2001-2004, we calculated population estimates using 
Zippin model M, (ZtPPiN, 1958; Мнть et al , 1982), which assumes a behavioral response to 
capture. Escaped salamanders were excluded from removal pass counts. Escapes were higher 
їп 1999 compared to 2001-2004, when the percent of salamanders that escaped per stream 
across all removal samples averaged 36 + 3.9 and 26 + 31°. for Desmognathus, 4% + 
73 sand 15 + 19 « for E. and 44 + 18 7 and 11 30 for Girmophilust 
Pseidotriton. respectively Basic assumptions for removal studies include population closure, 
equal sampling etfort, equal catchability. and effective reduction of the population after each 
search Unfortunately in the case of stream salamanders, some of these assumptions are 
difficult to meet (Вкуст, 1995) 


We used the “closed captures” selection m program MARK (Witt & BURNITAM, 1999) 
io calculate V., (А), p and s, (5) for capture-recapture model Mo. which assumes a constant 
capture probabil.ty, and removal model My, For Desmognathus data m 1999. we used paired 
1 tests to test for significant differences between the 7% of the capture-recapture and removal 
estimates and used program CONTRAST (Sauer & Wi Its. 1989: Hints & Sati R. 1990) 
to test whether рУ from the capture recapture and removal data diſſered among streams. For 
the 2001-2004 data, we compared whether fs estimated using removal models differed among 
streams within years and among years within streams for Desmognathi and Luriced using a 
Chi Square test implemented in program CONTRAST. Analyses were conducted within 
groups (streams or years), then pooled to provide a composite Chi Square test with summed 
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degrees of freedom. We also tested for differences among years and streams m À for 
Desmognathus and Eurycea using a two-way ANOVA in SPSS (Norusis, 1992) 


RESULTS 


Across the capture-recapture surveys at the six streams in 1999, we marked 180 Desmo- 
gnathus, 62 E bislineata, and 17 Gy rinoplulusi Pseudotriton (tab. I). Recapture rates were fairly 
low, ranging from 0 to 33 ? across species and stream sites (tab, 1), Because the numbers of 
marked and recaptured salamanders were low, the simplest model (model Mo), which 
assumes a constant capture probability, was usually the modcl of choice in program CAP- 
TURE; the results of this model are presented in tab 2. Unless a population 1s large or 
exhibits high capture probabilities, model selection may not be able to detect a pattern in js 
and will select the default model Mo (MENKENS & ANDERSON, 1988) Because capture- 
recapture estimates based on maximum-likelihood and Bayesian models were from the same 
data set, values for Ñ and р from these methods were quite similar; we used the estimates based 
on maximum-likelihood for analyses We were only able to calculate capture-recapture 
population estimates for E. bislimeata and Grrmophilusi Pseudotriton at one stream each, 
Jeremy's Run, where one individual of each species was recaptured (tab 2) Desmognathus 
individuals were recaptured at all 6 streams and capture-recapture fs averaged 0 06 + 0.014 
(range: 0 02 0 10) (tab. 2). Capture-recapture ÿ's differed among streams for Desmognathus 
( = 13.3, df — 5, Р = 0.02). Capture-recapture does not perform well unless is exceed 0.30 
(Ware et al., 1982), which was never the case, Because of this, the standard errors of & were 
sometimes very large (tab. 2). 


Based on the six pass removal data from 1999, population estimates could be calculated 
at all 6 stream transects for Desmognathus, 3 for E bishineata, and 2 for Gyrmophdus! 
Pseiudotriton (tab. 2). Removal ys averaged 0.25 + 0.077 (range: 0.08-0 61) for Desmognathus 
across stream transects, 0.25 + 0 079 (0 09-0 35) for E beat, and 0 44 + 0 060 (0.38- 
0 50) for GYrmoplulus Pseudotriton Removal p's differed among streams for Desmognathus 
(7 = 279, df S. P < 0.001) Removal js were significantly higher than those based on 
capture-recapture for Desmognathus {t = 2.2, df = 5, Р = 0.04, tab. 2). 


For the 2001-2004 data, we calculated removal population estimates for stream salaman- 
ders at stream transects at the eight streams (tab. 3), Estimation was not possible when zero 
counts occurred on the second pass when two passes were used or on the second and third 
passes when three passes were used, or when counts increased across subsequent passes, We 
could calculate population estimates at 55° , 54". and 13 «of the total 56 stream transects 
surveyed from 2001 to 2004 [or Desmognathus, Euryeca and Girmophilusi Pseudotriton, 
respectively (tub. 3). When estimates were calculable, p's averaged 066 + 0.025 (range 
0 39-0.96) for Devnognathie across stream transects. 0 64 + 0031 (0 27-0 89) for Eurycea. 
and 0 62 = 0065 (0 27 0 75) for GU Pscudotriton (tab. 3). These detection probabi- 
lues were much higher than those found in 1999 Lsing program CONTRAST, we found that 
P's diflered among streams for Pe litt (, — 29.7, df = 17, P – 0028) but not for 
Eurscea Detection probabilities did not duler among years at a stream for either Desmegne- 
thus or Eurycea 
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Table 1, - Numbers of salamanders summed across 5 to 6 capture-recapture surveys in 1999 that 
were too small to mark, escaped capture, or were marked or recaptured (% recaptured in 


parentheses). 
Species Stream ‘Not marked Escaped Marked Recoptured (%) 

Désmognatina Jeremy s Run E] са 5 BO) 
Keyser Run м 29 35 16) 
Land's Run E 30 2% 705) 
North Fork Thoenion 15 20 EI 20) 
Pass Run Tributary м 25 18 2010 
Piney River s 2 9 mm 

 Eurycea bishneata Jeremy's Run 7 14 3 Ime 
Keyser Run 4 ri £ 90 
Land's Run. 4 4 5 040) 
North Fork Thornton о 20 35 0 0% 
Pass Run Tributary 2 12 м 
Piney River y $ 3 00) 

Gyrinophilus / Pseudotriton | Jeremy's Run. 0 7 8 тазу 
‘Keyser Run ° t 1 00) 
Land's Run 0 3 1 LI] 
North Fork Thomion 0 1 2 00) 
Pass Run Tributary в ° 1 00) 
Pincy River 9 5 4 0 


Table 2. - Population estimates ( М + standard error, sx), 95 % confidence intervals (CI), detection 
probabilities (p+ з) and models used for species encountered during capture-recapture 
(CR) and removal (REM) sampling in 1999 at 6 streams m Shenandoah National Park. For 
Bayesian results, we present the mean № /mode №. 


Å (s) 95 %CI BG) | мма 
Desmognathus 

deremy's Run св в necs sers 010% % | o 
Bayesian | 6 124/119 259) 86185 | 0.100.023) 

REM 6 [mns 91004) 86131 | 025(0053) | в 

Land's Run св в 62(18.8) 4am | 009002) | о 
Bayssan | 6 72166 (25 6) 41-07 | 0090 030) 

REM 6 965854 62076) aus | 014(0.085) | в 

Pass Run Tributary CR E 71649) 31-202 | ооо | o 
Bayesso | 5 їз (1049) 32392 | 0060005) 

REM 6 1255659 108 (900) 52-537 | 0.08(0.083) | в 

Piney River CR 5 33094) 1366 | 0060030 О 
Bayesian | 5 | 95167 (1594) %% | 00000048 

REM 6 siaaa | 166 0230117) | в 

Keyser Run св 5 | 4934792) 120-2505 | 002004) | о 
Bayesian | 5 556419 461 4) 128-1844 | 0.020.015) 

REM 6 5431008 2100) 2-3 | овоот) | B 

North Fork Thomton CR 6 25301693) asn | 00000 | o 
Dyem | 6 316266 (275.5) 92104. | 003(00.6) 

REM 6 11.02.1746 70(18.3) sus | 0170972) в 

Eurycea butinoata 

Jeremy's Run CR 6 sen 320 | 01700139) | o 
Baysan | 6 19760 1) зао | 0120095) 

RLM в 50651 3) 2533 |owçns | в 

Pass Run REM 6 5a) 56 035016) | B 

Piney River. REM 6 199,10 300) 33 озб) | B 

Gyrimophitus Pseudotriton 

Jerry's Run cR в 27042) пыз? | оош) | o 
Bayesan | 6 } 24001384) 11433 [ 005006) 

North Fork Thornton км в | 100) 14 0500250 | в 

Piney River REM 6 1 300) эз osom |в 
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Table 3. Population estimates ( À + ). 95 % confidence intervals (CI) and detection probabilities 
(+ s) for salamanders encountered during removal samphng in Shenandoah National 
Park. Analyses were based on two or three diurnal passes conducted consecutively at one 
(2001) or two (2002-04) transects (T) at each of eight streams (escaped salamanders 
excluded) using model My (Zippin) from program CAPTURE. - indicates a third pass was 
not conducted No data for a particular species for a year. stream or transect indicates that 
none were detected or that estimation was not possible. 


Taxon - Stream Year Nis эз» Ci Pls) 
[Dewoanams 
Piney River zo 1 1038 22 0670272) 
Р ney fributany 2m fi 4030 44 06740 192) 
2 91 86) 920 0520245) 
з |! 30025) 33 07540217, 
он |! 302» 33 0750215 
2 340.27) 33 075 (0.217) 
Ivy Creek zor fi 44 (14 98) 34-109 0454020 
x2 |! 35079) 3535 090004 
2 4715.80) 4370 0630129) 
ю | 024389 1-34 0581020) 
2 2025 22 06710272) 
юм 1 106530) 934 039 (0215) 
2 90.09) 9s 0.69 (0 128) 
Doyle's River 201 |1 10(035) 10-10 091 (0087) 
2 2 124073) na 0 1 70 % 
20 1 (1.06) m 0.62 (0 135) 
Hawksbill жю |1 2801437 21-100 042 (0.278, 
2002 251021) 2525 0.90 % 038) 
2003 30) 33 06010219, 
2008 34071) 33 060 (0.219) 
40.51 з 073 (0 134) 
Pass Run 2002 9099 95 0040128 
2004 71087) M 064 (0 % 
70030 77 004 (01451 
Jeremy s Ran 2001 nom na 9 40 224 
2002 2004.15) 1840 046 (0 1773 
61038) 66 075(0153) 
200 1640 88) 1616 084 (00%, 
1866) 1351 051 (0307) 
200 160072) 16-16 0.73 00951 
10(0 63 1010 070102) | 
em 
Puy River x |2 von 08,0 % 
жз |2 20038) 0.06 (0272 
2w [1 non 04910215) 
Pane Tributan w (т am 06000193 
202 2 ET 031 0183) 
жю [1| 176. aom 0231014) 
2 | 127 2467) 054 (0228) 
me [ao 120073) 71011 
2] sat wam 0620121) 
[3973 sw [a | 71- 30020) 09 (0 105) 
Dav € $R ver w3 42. ed os 0750 159 
ws [i| en 1016 86) 00700 122 
[e эш fe) 6 ШЕЛ 0.740232) 
s | ава дио 0850007) 
zum |а| az- Gus 0750159 
ma t| n2 61 69) 0510294) 
2121 оа» ок уота 
tms Run aw 2444 ETEN O47 0 Mb, 
юы [1 | тозо ness O38 (01871 
2) 1497 жн бз a 
vm skan ma faje | sax) пот 
эю [ir] sio d th prm 
ma dius 07 EET 
2) 420 20 1 PETIRI 
unc à me Eau] изн 62447 27:628 
DTA uum ET ot 
2 ризо. reas n FEE 
icra рав чаки Ms Gor tse 
өн (ахо "mm о 0 5 40 l 
L 2j gres Sedem m [TENTE 
азон раль а cue! 
ue з | а2- e 120 «7500159 
p 1| en EE 19-10 EIS 
H m m nen, 
wes 1 mem EI cen эу 
зан HE МЕ 45 СЕИ 
а хт 19250 ma: OSL si 
2) 220 IET ЕЕ] 06719,9» 
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To analyze population change at a site, at least two years of data are needed and analyses 
should rely on population estimates to avoid bias associated with raw counts For Desmogna- 
thus and Eurycea, we had complete sets of for 2001-2004 at three streams each and À for 3 
of the 4 years at another two streams (tab. 3) We found significant differences in N across 
years (F = 12.7, df = 3.9, P = 0 001) and streams (F- 9 9, df = 4.9, P — 0 002) as well as a 
significant year*stream interaction (F = 4.6, df = 10,9, P = 0.015) for Desmognathus. but no 
significant differences for Eurycea Desmognathus population estimates were significantly 
higher in 2001 (24 + 7.8) and 2002 (20 + 5.4) compared to 2003 (9 = 2.7) and 200448 + 1 4) 
The park experienced heavy precipitation during the summer of 2003, with average stream 
flow rates 19 and 7 times higher than stream flow rates in the summers of 2002 and 2001. 
respectively (Shenandoah Watershed Study data, Rick Webb, pers comm.) Data from July 
2004 are not yet available, but 2004 flow rates were most likely termediate between the 2003 
and 2001-2002 flow rates. 


DISCUSSION 


Long-term monitoring programs require cost-effective and efficient techniques to gather 
accurate and precise data. Unfortunately, the spatially variable (i. e., significant differences 
among streams) and sometimes low detection probabilities found in this study using capture- 
recapture and removal methods reinforce the need for estimating 5's as part of stream 
salamander abundance estimation studies. Our study also indicates the importance of deve- 
loping better methods for estimating stream salamander populations such that estimates are 
consistently available on a yearly basis for trend analyses. 


We found that removal sampling yielded higher p's for stream salamanders than capture- 
recapture sampling. Other capture-recapture surveys of stream salamanders have also shown 
low recapture rates and hence detection probabilities (BaRTHALMLS & Bi (us, 1972, Nun ts 
& KAPLAN, 1998). Indeed. Макоо (2001) used VIE to mark 44 Е bisimeata and D fuscus 
but did not recapture any in the field. Bruce (1995) used removal sampling (7 passes set 2-3 
days apart) and found low to moderate standard errors for population esumates of D 
monticola and suggested removal sampling was a promising technique to monitor salamander 
demographics. Other factors favoring removal over capture recapture sampling are that 
removal sampling usually requires shorter sampling intervals, reduced field personnel, 
and less funding than capture-recapture, and appears to be ideal for amphibians such as 
aquatic larvae that are highly detectable and have limited home ranges and mobility (Наук. 
1994). 


If removal sampling 15 to be used for long-term monitoring. field protocols play an 
important role in determining their success. In our removal surveys, effective reduction of 
populations sometimes did not occur even after six passes. This may be due m part to the high 
percentage of salamanders that escaped capture, though if we analyzed removal data includ- 
ing escapes in the passes, the same issues would be apparent. It is important to note that when 
the percent of escapes were lower as they were in 2001 2004 compared to 1999, the p's for 
species were higher. With fewer escapes and larger sample sizes. there 1s potential for better 
estimates. Removal estimates using Zippin's method are unreliable if less than half the 
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population is removed (Bruce, 1995), but WHITE et al. (1982) considered detection probabil- 
ities greater than 0.20 adequate for estimating population abundance in removal experiments. 
BRucr (1995) found that 7 passes probably reduced total D monticola populations by more 
than half at his study sites, but he had difficulties reducing numbers of first year Juveniles, 
which may have shown “increased surface activity ..as the larger salamanders were removed 
(i.e, a response to reduced competition or predation)”. SOL THERLAND et al. (2004) used 
two-pass removal sampling and were unable to calculate population estimates for species at an 
average of 75 % of the streams surveyed because salamander numbers did not decrease or 
were zero in the second pass. Removing salamanders from under the top layer of rocks may 
disturb or “unearth” other salamanders deeper in the rock substrate As we sometimes 
observed, this can lead to more salamanders in the surface population during subsequent 
passes than in the first pass before disturbance. 


Several factors could be changed in our removal protocol to improve Ñ and f estimates, 
Conducting surveys on wet or humid nights, when more of the salamander population may be 
on the surface foraging, might yield better removal estimates. D fuscus and Е bistineata 
emerge one hour after sunset (HOLOMUZK1, 1980) and D monticola emerge shortly after dark, 
with peak activity occurring around midnight and again at dawn (SHEALY, 1975; HAIRSTON, 
1986). However, working at night along rocky streams can be difficult and treacherous. 
Another option would be to conduct more removal passes, providing the option to group data 
from earlier passes in which no decreases in removals occurred Pilot studies m which a large 
number of removal passes are conducted to determine the appropriate number and grouping 
of passes may be useful Another factor to consider is the size and placement of transects or 
plots. In our surveys, we only searched narrow 1- or 2-m bands along and/or in the stream. 
Most stream salamanders move between the stream channel, splash zone and bank Home 
ranges of D fuscus have been shown to vary tremendously, from 1.4 m? in Ohio (ASHTON, 
1975) to 25-114 m° in Kentucky (BARBOUR et al. 1969) D. monticola home ranges were 
estimated to be 8 4 m? m Kentucky (HARDIN et al., 1969). During warm months, E hishinearu 
tagged with radioactive isotopes moved within a 14 m° area (ASHTON & Аѕнтом, 1978), but in 
June some post-breeding migrants moved more than 100 m from a stream (MACCULLOCH & 
Biprr, 1975), which probably explains the particularly low recaptures we observed for this 
species m the 1999 capture-recapture surveys. Surveying a wider area of bank along with the 
stream channel to incorporate more of the target species’ individual home ranges may yield 
better removal estimates. 


Other new approaches may prove to be more useful for stream salamander population 
estimation, Our removal estimates were based on populations at single stream transects. New 
analytical methods developed by Rovt.r (20044-^). Royir et al (2001) and Dorazio et al (in 
press) aggregate information across sample sites such that removal sampling can estimate the 
abundance of spatially distinct subpopulations These models incorporate spatial models of 
abundance (eg. Poisson, negative binomial) with models of detection probability and have 
been shown to yield abundance estimates with “similar or better precision than those 
computed using the conventional approach of analyzing the removal counts of each subpop- 
ulation separately" (Dorazio et al., in press). 


A different approach would be to estimate the proportion of area (in this case, streams) 
occupied (PAO) by stream salamanders over time (MACKENZII et al., 2002), The РАО method 
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estimates site occupancy and detectability of species based on presence/absence data recorded 
from repeated visits to sites selected using a probabilistic sampling frame within an area of 
inference. Stream salamander species that exhibit low detection probabilities and occupy 
fewer sites would require more streams and visits per stream for PAO estimation (MACKENZIE 
& RovLE, submitted). Note that repeated visits to streams could be satisfied by surveying 
multiple transects along the length of a stream. 


Despite the problems evident in this study, population estimation efforts incorporating 
detection probabilities may be necessary to assess trends in stream salamander populations, 
Better survey methods (e.g., transect designs) and population estimation techniques (e.g... 
aggregated removal or PAO approaches) need to be tested and developed such that reasonably 
low bias population estimates can be consistently calculated for sites over ume. In addition, 
spatial design of sampling associated with hypothesis testing incorporating covariates that 
may influence stream salamanders (р, N, site occupancy), such as the percent of impervious 
surface in a watershed and stream flow rates, should be incorporated alongside monitoring to 
best yield inferences about how changes in stream salamander populations over time are 
influenced by environmental factors. 


RÉSUMÉ 


Dans l'est des États-Unis, les salamandres torrenticoles de la famille des Plethodontidae 
représentent une biomasse élevée dans et auprès des ruisseaux issus des sources, Elles peuvent 
ainsi constituer d'intéressants indicateurs de la santé de ces écosystèmes. Beaucoup d'etudes 
de ces salamandres se sont appuyces sur des indices démographiques utihsant des décomptes 
d'animaux et non pas sur des estimations fondées sur des techniques comme les captures- 
recaptures ou le ramassage des individus. L'emploi de procédures d’evaiuation permet le 
calcul de probabilités de détection (la proportion d'animaux réellement presents détectés lors 
d'uneetude)et de leur ecart-type, et peut permettre de déterminer les tailles et les dynamiques. 
des populations de salamandres. En 1999, nous avons employé les méthodes de capture- 
recapture et de ramassage pour évaluer des populations de salamandres du genre Desmogna- 
thus dans six ruisseaux du Parc National de Shenandoah (Virginie, Etats-Unis) La méthode 
du ramassage s'est avérée plus efficace elle a donné des probabilites de détection plus élevees 
que celle de capture-recapture Pendant la période 2001-2004. nous avons employé la methode 
du ramassage dans et auprès de huit ruisseaux du Pare afin d'évaluer la fiabilite de cette 
technique pour la surveillance à long terme de ces populations de salamandres Lors de 
transects le long des ruisseaux, nous avons obtenu des probabilites de detection de 0.39 à 0.96 
pour Desmognathus, de 0.27 à 0 89 pour Eurycea et de 0,27 à 0,75 pour Gyrmophilus 
porphyrincus/Pscudotriton тибет Les probabilités de détection n'ont pas varié au cours des 
annees pour Desmognathus et Euricea, mais ont différé selon les ruisseaux pour Desmognd- 
thus Les evaluations des populations de Desmognathus ont diminue entre 2001-2002 et 
2003-2004, ce qui peut être hé à des changements dans le régime hydrique des ruisseaux. Les 
procedures de ramassage constituent une methode fiable pour l'évaluation de populations de 
ces salamandres, mais les méthodes de terran doivent être conçues de mamère à remplir les 
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conditions statistiques des méthodes d'échantillonnage. De nouvelles méthodes d'estimation 
des populations de ces salamandres sont discutées. 
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